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ABSTRACT 



Ultra-luminous Compact X-ray Sources (ULXs) in nearby spiral galaxies and Galac- 
tic superluminal jet sources share the common spectral characteristic that they have 
unusually high disk temperatures which cannot be explained in the framework of the 
standard optically thick accretion disk in the Schwarzschild metric. On the other hand, 
the standard accretion disk around the Kerr black hole might explain the observed high 
disk temperature, as the inner radius of the Kerr disk gets smaller and the disk tem- 
perature can be consequently higher. However, we point out that the observable Kerr 
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disk spectra becomes significantly harder than Schwarzschild disk spectra only when 
the disk is highly inclined. This is because the emission from the innermost part of 
the accretion disk is Doppler-boosted for an edge-on Kerr disk, while hardly seen for 
a face-on disk. The Galactic superluminal jet sources are known to be highly inclined 
systems, thus their energy spectra may be explained with the standard Kerr disk with 
known black hole masses. For ULXs, on the other hand, the standard Kerr disk model 
seems implausible, since it is highly unlikely that their accretion disks are preferen- 
tially inclined, and, if edge-on Kerr disk model is applied, the black hole mass becomes 
unreasonably large (> 300Mq). Instead, the slim disk (advection dominated optically 
thick disk) model is likely to explain the observed super-Eddington luminosities, hard 
energy spectra, and spectral variations of ULXs. We suggest that ULXs are accreting 
black holes with a few tens of solar mass, which is not unexpected from the standard 
stellar evolution scenario, and that their X-ray emission is from the slim disk shining 
at super-Eddington luminosities. 

Subject headings: superluminal jet sources: ultra-luminous X-ray sources: accretion 
disks, slim disks: Schwarzschild black holes, Kerr black holes 

1. Introduction 

Ultra-luminous compact X-ray sources (ULXs) have been found in nearby spiral Galaxies, with 
typical 0.5 - 10 keV luminosities 10 39 to 10 40 erg s _1 (e.g., Fabbiano 1988; Colbert and Mushotzky 
1999; Makishima et al. 2000; Colbert and Ptak 2002; Foschini et al. 2002). These luminosities are 
too small for AGNs, and most ULXs are in fact located significantly far from the photometric center 
of the galaxies. On the other hand, ULXs are too luminous to be considered as the same class of 
the compact binary X-ray sources in our Galaxy whose luminosities are almost always < 10 39 erg 
s" 1 . 

Significant time variations have been detected from ULXs, and their energy spectra are suc- 
cessfully modeled with emission from optically thick accretion disks (Okada et al. 1998; Mizuno et 
al. 1999; Kotoku et al. 2000; Makishima et al. 2000; Mizuno, Kubota and Makishima 2001), as is 
the case for the "High" state (= Soft-state) of Galactic black hole candidates (e.g., Tanaka and 
Lewin 1995). In addition, discovery of bimodal-type spectral transitions (Kubota et al. 2001) and 
orbital modulations (Bauer et al. 2001; Sugiho et al. 2001) from several ULXs further demonstrate 
their resemblance with Galactic black hole candidates. These observational facts suggest that ULXs 
are moderately massive black holes, which may be scale-up versions of the Galactic black holes. 
So that the observed luminosities of ULXs, ~ 10 40 erg s _1 , do not exceed the Eddington limit 
Lem = 1-5 x 10 38 (M/Mq) erg s -1 , the black hole mass has to be as large as or greater than ~100 
Mq, assuming isotropic emission. How to create such "intermediate" mass black holes, if truly 
exist, is an intriguing question (e.g., Ebisuzaki et al. 2001). On the other hand, if X-ray emission is 
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unisotropic, the black hole mass is not required to be so large. For example, a beaming model for 
ULXs has been proposed (King et al. 2001; Kording, Falcke and Markoff 2002); however, presence 
of a bright nebula surrounding M81 X-9 (Wang 2002) suggests that the nebula is powered by the 
central X-ray source and that strong X-ray beaming is rather unlikely. Another possibility to solve 
the super-Eddington problem through unisotropic emission is a geometrically thick accretion disk 
(Watarai, Mizuno and Mineshige 2001), which we favorably consider later (section 4.1). 

GRS 1915+105 and GRO J1655-40 are the two well-known Galactic superluminal jet sources 
and established black hole binaries with reliable mass measurements (M6.3 ± 0.5M Q for GRO 
J1655-40 [Green, Bailyn and Orosz 2003] and « 14M for GRS 1915+40 [Greiner, Cuby and 
McCaughrean 2001]). If soft X-ray energy spectra of these sources are fitted with an optically thick 
accretion disk model, their characteristic disk temperatures are ~1.3 - 2.0 keV (e.g., Belloni et 
al. 1997; Zhang et al. 1997; Zhang, Cui and Chen 1997; Tomsick et al. 1999). These values are 
systematically higher than those of ordinary and well-studied soft-state black hole candidates such 
as Cyg X-l and LMC X-3, whose disk temperatures are almost always less than ~ 1 keV (e.g., 
Tanaka and Lewin 1995). 

Okada et al. (1998) found that a luminous ULX in IC342 (~ 2 x 10 40 erg s _1 ; "Source 1") 
has unusually high disk temperature (~1.7 keV), which is similar to the accretion disk spectra of 
Galactic superluminal jet sources, and pointed out that such a high disk temperature cannot be 
explained in the framework of the "standard" accretion disk around a Schwarzschild black hole 
(see section 2.1 and 3.2). The "standard" disk denotes the situation that all the gravitational 
energy release is converted to thermal radiation, in contrast to the "slim" disk which is an optically 
and geometrically thick disk with dominant energy advection (section 4.1). Makishima et al. 
(2000) summarized ASCA observations of seven ULXs, and concluded that the unusually high 
accretion disk temperature is a common spectral property of ULXs. King and Puchnarewicz (2003) 
pointed out that ULXs, as well as ultrasoft AGNs, violate the apparent blackbody-temperature and 
luminosity relationship which is required not to exceed the Eddington limit. Zhang, Cui and Chen 
(1997) and Makishima et al. (2000) suggested that the unusually high disk temperature of the 
Galactic superluminal jet sources and ULXs may be explained if they harbor fast rotating black 
holes, since inner edge of the accretion disk gets closer to the black hole in the Kerr geometry 
and the disk can be hotter. On the other hand, Watarai et al. (2000) and Watarai, Mizuno and 
Mineshige (2001) proposed that, instead of the standard disk, the slim disk model may explain the 
X-ray energy spectra of Galactic superluminal jet sources and ULXs. 

In this paper, we focus on the "too-hot disk" problem of ULXs and Galactic superluminal 
jet sources. We apply standard accretion disk models in Newtonian, Schwarzschild, and Kerr 
cases, and discuss how the relativistic effects and the disk inclination affect the black hole mass 
and the mass accretion rates obtained from the model fitting. We see that observed hard spectra 
of Galactic superluminal jet sources may be explained by strong relativistic effects in the highly 
inclined standard Kerr accretion disk. On the other hand, the hard spectra and the super-Eddington 
problem of ULXs are difficult to explain in the framework of the standard accretion disk model. 
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We see that the slim disk model may explain the observed super-Eddington luminosities and hard 
energy spectra of ULXs more naturally 

2. Characteristics of the Standard Accretion Disk Model 

Before quantitative discussion of the observed energy spectra and accretion disk parameters, 
we summarize important observational characteristics of the standard accretion disk model. 

2.1. Disk Temperature 

Let's assume a standard optically thick accretion disk (Shakura and Sunyaev 1973) in the 
Schwarzschild metric, in which case the last stable orbit around the non-rotating black hole is 6 r g , 
where the gravitational radius r g = GM/c 2 is defined. We identify the last stable orbit with the 
innermost disk radius, rj n . The radial dependence of the "effective" temperature of an optically 
thick accretion disk may be written as, 

m , , f 3GMM n . , . 1 1/4 

Teff(r) = ^^ r R R (r/r m )^ , (1) 

where Rr is ^1 — y/ri n /rj in the Newtonian case, otherwise includes additional general relativistic 
correction (see e.g., Krolik 1999). The effective temperature is zero at r = rj n , peaks at r « 8 r g , 
and decreases with oc r~ 3 / 4 outward. 

Because of the Comptonization in the disk atmosphere, the color temperature of the local 
emission becomes higher than the effective temperature, and their ratio is almost constant at ~ 1.7 
(section 2.3). Therefore, the observable disk color temperature also peaks at r « 8r g , and takes 
the maximum values, 

e -»,,o-,3,w(^M)(*) 1/4 (^)- 1/4 , , 2 ) 

where we define the critical mass accretion rate M c = 2.9 x 1O 18 (M/M ) g s _1 , so that M = M c 
gives the Eddington luminosity (see appendix A). The small temperature variation reflects the 
inclination effects (from face-on to i = 80°), such that inclined Schwarzschild disks have slightly 
harder spectra (section 2.4 and Appendix B). 

From equation (2), we can see that an optically thick accretion disk around a 7 M & or 100 
Mq Schwarzschild black hole may not have a higher color temperature than ~1.3 keV or ~0.7 keV, 
respectively, unless super-Eddington mass accretion is allowed. To explain higher disk temperatures, 
either super-Eddington mass accretion rates or unreasonably small mass is required. This is the 
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"too-hot disk" problem of the Galactic superluminal jet sources and ULXs we are concerned in this 
paper (see also King and Puchnarewicz 2003). 

2.2. Multi-color disk approximation 

We point out that the maximum disk color temperature as expressed with (2) is a quantity 
directly constrained from observations. If we approximate optically thick accretion disk spectra 
with a simple multicolor disk blackbody model (MCD model; Mitsuda et al. 1994; Makishima et al. 
1996) in which local emission is a blackbody and radial dependence of the temperature is simplified 
as T(r) = Ti n (Rin/r) 3 ^ , the apparent inner disk radius R{ n and temperature T« n will be the 
independent model parameters 2 , such that the spectral shape is determined only by Tj n , and the 
flux is proportional to R 2 n Tf n . When observed spectra are fitted with the MCD model, Tj„ may be 
identified with the maximum disk color temperature (2), while the true inner disk radius ri n and 
apparent radius Ri n have a relation ri n rs 0.4 (T co i/T e ff) 2 Ri n (Kubota et al. 1998). 

2.3. Local Energy Spectra 

In the inner part of the optically thick accretion disk, the electron scattering opacity is dom- 
inant, and the hot disk atmosphere distorts the emergent spectra via Comptonization. Precise 
theoretical calculation of the accretion disk spectra for Galactic black hole candidates has been 
made by several authors (e.g., Shimura and Takahara 1995; Ross and Fabian 1996; Blaes et al. 
2001). These authors tend to agree that the local spectrum from each ring of the disk can be ap- 
proximated by the "diluted blackbody", (T e ff/T co i) 4 B(E,T co i), in the ~0.5 - 10 keV band, where 
B(E,T co i) is the Planck function with color temperature T co \ (> T e ff). Also, ratio of the color 
temperature to the effective temperature, T co i/T e ff, is virtually constant at 1.7 - 1.9 along the disk 
radius for high values of the accretion rate and low values of the viscosity parameter (Shimura and 
Takahara 1995). Hence, as long as T co i/T e ff is assumed constant, equation (2) is valid and simple 
MCD approximation may be used to describe observed disk spectra. Throughout this paper, we 
assume T co i/T e ff = 1.7 and the diluted blackbody local spectra (see discussion in section 3.3.1 for 
the limitation of this assumption). 

It should be noted that the mass obtained by fitting the standard accretion disk model to the 
observed spectra is proportional to (T co i/T e ff) 2 , as long as the diluted blackbody approximation is 



In this p&pGr, Tin denotes the real disk inner radius, while Ri n is the apparent inner radius as a MCD model 
parameter. Appendix C gives conversion formulae between the MCD parameters (Ri„ and T„) and realistic disk 
parameters (M and M) in the Schwarzschild case. 
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valid. This can be readily seen as follows: The accretion disk spectrum may be expressed as, 

2 ^ r B ( E i T col(r,M,M)} dr 

~r in 2 !°° 2i:{r/r in )(^A B (E,T col (r/r m , M, M)) d(r/r in ) 

J{r/r in )=l \ J- col J ^ ' 

= f 6G ^ff/^oi) 2 M \ r°° 2 7r(r/r in )B(E,T eff (r/r in ,(T eff /T ml ) 2 M,M)) d(r/r in ), (3) 

V c / J(r/r in )=l V > 

where we have used rj n = 6GM/c 2 and T co \ oc M 1//4 (r/rj„)~ 3 ^ 4 ((T e jj/T co ;) 2 M) from equation 
(1). Equation (3) indicates that the accretion disk spectrum with the black hole mass M and 
the diluted blackbody local spectrum is identical to the disk spectrum with the black hole mass 
(T e ff /T co i) 2 M and the blackbody local spectrum. 



2.4. Relativistic Effects 

Due to the relativistic effects in the Schwarzschild metric, observed disk spectral shape is mildly 
inclination-angle dependent, such that inclined disk spectra become slightly harder (appendix C, 
figure 1). This is due to the gravitational redshift and Doppler boosts near the inner edge of the 
accretion disk. Note that the relativistic effect is not very significant outside of 6r 9 , such that the 
Schwarzschild disk spectra are not very much different from those of Newtonian disks. This is in 
contrast to the extreme Kerr case, in which the relativistic effect is enormous necir t^yi ~ r a and 
the disk spectrum is extremely inclination dependent (see below; figure 1 and 2). 

While in the Schwarzschild metric the last stable orbit around the black hole is 6 r g , in the Kerr 
metric it can go down to f .24 r g with an extreme angular momentum of a = 0.998. As the inner 
edge of the accretion disk approaches the black hole, more gravitational energy is released, and the 
inner disk temperature can get higher. X-ray energy spectra from optically thick Kerr accretion 
disks have been calculated by many authors including Cunningham (1975), Connors, Piran and 
Stark (1980), Asaoka (1989), Sun and Malkan (1989), Laor, Netzer and Piran (1990), Hubeny et 
al. (2000) and Gielihski, Maciolek-Niedzwiecki and Ebisawa (2001). In this paper, we calculate the 
Kerr disk spectral model using the transfer function by Laor, Netzer and Piran (1990) for a = 0.998. 
To make a comparison easier, we assume rj n = 1.24 r g , and the same local spectrum as we did for 
the Schwarzschild disk model (section 2.3), namely diluted blackbody with T co i/T e ff = 1.7. 

In Figure 1, we show comparison of the Schwarzschild and Kerr disk spectra, as well as New- 
tonian one, for the face-on and a highly inclined case (i = 80°). The biggest difference between the 
Kerr disk model and the other two models is the presence or absence of the innermost region of the 
disk from 1.24 r g to 6r g . When the disk is close to face-on, the contribution from this part is not 
significant because of the gravitational red-shift and light-bending; as a result, the total Kerr disk 
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spectrum is not very much different from the Newtonian or Schwarzschild ones. On the other hand, 
emission from the innermost part becomes very significant for a highly inclined disk due to Doppler 
boosting. Consequently, a near edge-on Kerr disk spectrum comprises a significant amount of high 
energy X-ray photons compared to the Newtonian or Schwarzschild disks. 

Inclination angle dependence of the observed Kerr disk spectrum is different in different energy 
bands (figure 2). At lower energies, where most of the emission is from the outer parts of the disk 
(r > 400 r g ), the flux is proportional to the projected disk area and decreases with inclination, just 
as Newtonian disks. On the other hand, higher energy fluxes from the innermost parts of the disk 
(1.26 r g < r < 7r g ) are enhanced as the disk is more inclined, due to strong Doppler boosting. 
Consequently, inclined Kerr disks are brighter than the face-on disks in the highest energy bands, 
as previous authors have already pointed out. 

3. Application of the Standard Disk Model 

We apply the standard accretion disk spectral models to ASCA archival data of GRO J1655-40 
and IC342 Source 1. The observations were made on 1995 August for GRO J1655-40 (when the 
source is in a bright state) and 1993 September for IC342 Source 1 (when the source is in the "high" 
state; Kubota et al. 2001; Kubota, Done and Makishima 2003). These datasets may be considered 
"exemplary", so that the same spectral data of GRO J1655-40 have been analyzed in Zhang et al. 
(1997), Zhang et al. (2000) and Gielihski, Maciolek-Niedzwiecki and Ebisawa, K. (2001), and those 
of IC 342 Sourcel in Okada et al. (1998), Mizuno et al. (1999), Watarai, Mizuno and Mineshige 
(2001) and Kubota et al. (2001). 

We assume the diluted blackbody local spectra with T co [/T e ff = 1.7 in sections 3.1 and 3.2. 
Difference among the Newtonian, Schwarzschild, and Kerr disk models are only the relativistic 
effects, including the variation of the innermost radius and energy conversion efficiency (appendix 
A). The Schwarzschild disk model we use (GRAD model) is explained in Hanawa (1989) and 
Ebisawa, Mitsuda and Hanawa (1991) (see also appendix B). The GRAD model is available in the 
XSPEC spectral fitting package, and this model also provides the Newtonian disk spectra we use 
in this paper, when the relativistic switch is turned off. 

3.1. Application to GRO J1655 40 

In Table 1, we show results of spectral fitting of the ASCA GIS spectra with Newtonian, 
Schwarzschild, and extreme Kerr disk models. System parameters have been determined from 
optical observations (Green, Bailyn and Orosz 2003), such that distance is 3.2 kpc, inclination 
angle is 70°, and M = 7M . We fixed the distance, while several different inclination angles are 
tried from 0° to 80°. The energy spectrum has a power-law hard-tail but contribution of which is 
minor below 10 keV and hardly affects discussion of the accretion disk spectrum. The power-law 
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slope is fixed at 2.5 which matches the simultaneous high energy observation with BATSE (Zhang 
et al. 1997), and its normalization is allowed to be free within the range acceptable by BATSE. 
Hence, the free parameters are M, M, power-law normalization, and the hydrogen column density. 
Besides minor local features (see Gielihski, Maciolek-Niedzwiecki and Ebisawa 2001), all the models 
with different inclination angles can reasonably fit the observed accretion disk spectrum. 

Note that the Schwarzschild disk model at the correct inclination angle (i = 70°) gives a too 
small mass (1.8M ) which is unacceptable (red curve in Figure 3 left). If we fix the mass at the 
correct value (7M ), the disk spectrum cannot be hard enough to explain the observed spectrum 
(the "best-fit" is shown with green curve in Figure 3 left): The "too- hot" disk problem is thus 
evident (equation 2). 

It is of interest to see how the black hole mass depends on the inclination angle. Observation can 
constrain the projected disk area, which is proportional to square of the mass. Hence M oc 1 jyj cos i 
in the Newtonian case, and the black hole mass at i = 80° is 2.4 times higher than that for the 
face-on disk. On the other hand, the mass increasing factor is 2.9 in the Schwarzschild case, 
and as much as ~ 10 in the extreme Kerr case. This is because inclined relativistic disk spectra 
become harder, which is most conspicuous in the extreme Kerr disk (section 2.4; figure 1). Given 
the observed spectra, to compensate the spectral hardening with inclination, the mass becomes 
necessarily larger (equation 2). 

The extreme (a=0.998) Kerr disk model with the correct inclination angle gives M = 15.9M . 
Compared to the fit with Schwarzschild disk model (1.8M ), significant increase of the mass does 
indicate the spectral hardening of the Kerr disk model. In fact, the derived mass is too large to be 
consistent with the realistic mass 7 M , which probably suggests that a=0.998 is too high. In fact, 
Gielihski, Maciolek-Niedzwiecki and Ebisawa (2001) applied the Kerr disk model to the same ASCA 
energy spectrum of GRO J1655-40, and concluded that a is likely to be between 0.68 and 0.88 to 
be consistent with M = 7M . To summarize, in agreement with Gielihski, Maciolek-Niedzwiecki 
and Ebisawa (2001), a standard Kerr disk model may explain the observed accretion disk spectra 
of GRO J1655-40 at the inclination angle i = 70°. The Kerr metric is required, but not with the 
maximum angular momentum. 

That GRO J1655-40 has probably a spinning black hole is also suggested by the recent dis- 
covery of 450 Hz QPO (Strohmayer 2001), as this frequency is higher than the Kepler frequency 
at the innermost stable orbit of a 7M Schwarzschild black hole. Precise analysis of the QPO 
characteristics from GRO J1655-40 suggests that the black hole is in neither a Schwarzschild nor 
a maximal Kerr black hole (Abramowicz and Kluzniak 2001), which is consistent with our energy 
spectral analysis. 
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3.2. Application to IC342 Source 1 

We assume the distance to the source 4 Mpc (Okada et al. 1998 and references there in). 
If isotropic emission is assumed, the luminosity will be 1.7 x 10 40 erg s _1 (1 - 10 keV) at this 
distance, thus M > 100M Q is expected so as not to exceed the Eddington luminosity. In table 2, we 
summarize results of the spectral fitting of the ASCA GIS spectrum with Newtonian, Schwarzschild, 
and Kerr disk models, for different inclination angles. All the fits are acceptable, thus models may 
not be discriminated based on the quality of the fits. Note that the disk luminosities always exceed 
the Eddington limit more than 10 times in the Schwarzschild case. 

The right panel in figure 3 indicates the ASCA GIS spectrum and the best-fit Schwarzschild 
disk model (in red). For comparison, we show a Schwarzschild disk spectrum with M = 1OOM 
at the Eddington limit (in green); we can clearly see that such an accretion disk has a too low 
temperature to explain the observed hard spectrum. 

The face-on Kerr disk model gives M = 27.3M and M = 2.0 x 10 20 g s" 1 (M = 16 M c ). A 
factor of ~ 3 increase of the mass compared to the Schwarzschild case is due to slight hardening of 
the face-on Kerr disk spectrum. Note that the super- Eddington problem does not disappear, as the 
face-on Kerr disk spectrum is not very different from the Schwarzschild one (section 2.4; figure 1). 
If we assume a very inclined disk with i = 80°, we obtain M = 332M and M = 1.5 x 10 20 g s _1 
(M = 1.0 Mc)- Now the super-Eddington problem is solved, that is a consequence of the fact that 
the inclined Kerr disk spectrum is much harder than the Schwarzschild one. However, there will 
be two serious problems to accept the inclined Kerr disk model for ULXs in general: First, we do 
not know a mechanism to create ~ 35OM black holes. Second, it is very unlikely that most of the 
accretion disks in ULXs are largely inclined when seen from the earth. 

3.3. Examination of Variants of Standard Disk Model 

In this section, we critically examine two possibilities which have been proposed to make the 
standard accretion disk spectra look harder and to explain the observed hard spectra of ULXs and 
Galactic super-luminal jet sources. 

3.3.1. Large T col /T eff ? 

A naive solution to explain the apparently hard disk spectrum may be to allow T co i/T e ff to be 
much greater than the standard value ~ 1.7. In this case the color temperature of the disk can be 
much higher for the same mass and mass accretion rate, while the effective temperature remains the 
same (equation 2). In other words, for a given observed disk flux and spectrum, M oc {T co i/T e ff) 2 , 
thus a larger mass is allowed (equation 3). 
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In fact, in the case of GRO J1655-40, (Section 3.1, table 1), if T col /T eff = 1.7y/7M Q /1.8M Q = 
3.4, the observed spectrum is explained with the Schwarzschild disk model with 7M & for the 
same distance and inclination angle. Borozdin et al. (1998) reached a similar conclusion that 
Tcoi/T e ff = 2.6 is required for GRO J1655-40 to fit with a Newtonian disk model whose inner 
disk radius is 3 times the Schwarzschild radius. For IC342 Source 1, we need a still higher value 
of T co i/T e ff = 1.7y / lOOM Q /8.9M0 =5.7, to fit with a face-on Schwarzschild disk model around a 
1OOM black hole. 

In this manner, extremely high T co i/T e ff values might solve the "too-hot" disk problem. In fact, 
it is suggested that standard disks with high viscosity may require such high values of the spectral 
hardening factor (Shimura and Takahara 1995). Given our lack of firm theoretical understanding 
of the mechanisms of viscous angular momentum transport (e.g., Merloni 2003), it would be better 
not to completely exclude the possibility of extremely large values of T co i/T e ff at this moment. 

On the other hand, there are several circumstantial evidences that T co i/T e ff cannot signifi- 
cantly exceed ~ 2: Independent theoretical calculations agree with values of T co i/T e f f = 1.7 - 1.9 
for accretion disks around ~ lOM black hole for high values of the accretion rate and low values 
of the viscosity parameter (Shimura and Takahara 1995; Ross and Fabian 1996; Blaes et al. 2001). 
Values of T co i/T e jf increase very slowly with mass, but cannot be higher than ~ 2.5 even for a 
1O 6 M black hole (Ross, Fabian and Mineshige 1992; Blaes et al. 2001). For observed accretion 
disk spectra of Galactic black hole candidates or weakly magnetized neutron stars, T co i/T e ff = 1.7 
—1.9 gives reasonable mass values which are consistent with those determined from dynamical mea- 
surements (e.g., Ebisawa, Mitsuda and Hanawa 1991; Ebisawa et al. 1993; Shimura and Takahara 
1995; Dotani et al. 1997). 



3.3.2. Comptonization of disk photon? 

If a standard optically thick disk is shrouded by a thermal corona whose temperature is slightly 
higher than the disk temperature, significant parts of the disk photons are comptonized and appear 
in higher energy bands. Thus, in such a situation, the observed hard energy spectra from Galactic 
superluminal jet sources and ULXs might be explained. 

The comptonization model may phenomenologically fit the observed spectra of GRO J 1655-40 
(Zhang et al. 2000; Zycki, Done and Smith 2001). Fitting with a thermal comptonization model (in 
which seed photons have a single temperature blackbody spectrum) gives typically the blackbody 
temperate < 0.5 keV, plasma temperature ~ 1.5 keV and scattering optical depth ~10. At such 
a low temperature and high optical depth, the heavy elements may not be fully ionized, and the 
plasma can be optically thick to photoelectric absorption. Presumably, such a physical situation is 
more reasonably described by standard optically thick accretion disk in which emergent spectrum 
is distorted by Comptonization (T co i > T e ff), which successfully explain the observation (section 
3.1). 
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We applied the comptonization model to IC342 Source 1 in 1993. When the seed photon 
is assumed to be disk blackbody, we obtain Tj n ~ 1.2 keV, plasma temperature ~ 2.5 keV, and 
scattering optical depth ~ 10 (x 2 /dof=51.8/46). Note that the Tj n cannot go down to ~ 0.7 keV as 
desired for the Schwarzschild disk around a 100 M black hole. Therefore, the original "too-hot" 
disk problem may not be solved even introducing the ad hoc comptonization plasma to harden 
the standard disk spectrum. When the seed photon spectrum is assumed to be blackbody, the 
blackbody temperature is 0.52 keV, the plasma temperature is 1.8 keV and the optical depth is ~ 
18 (x 2 /dof=51.8/46). Although this model may fit the data, it is hard to interpret these parameters 
in the physical context, and the super-Eddington problem is unanswered. Presumably, the slim 
disk model explained in the next section is more likely for ULXs when they are in the bright state. 
Significant disk comptonization may be taking place, instead, when ULXs are much dimmer and 
their energy spectra are power-law like (Kubota, Done and Makishima 2003). 



4. The Slim Disk Model 

4.1. Characteristics of the slim disk model 

So far, we have considered the standard optically thick accretion disk (Shakura and Sunyaev 
1973) in which radial energy advection is neglected and all the gravitational energy release is 
converted to thermal radiation. In accretion disk theory (for a review, e.g., Kato, Fukue and 
Mineshige 1998), there is another stable optically thick solution, which is called "optically thick 
ADAF (advection-dominated accretion flow)" disk or, "slim" disk (Abramowicz et al. 1988), which 
takes place when mass accretion rate is around the super-Eddington rate or higher. In the slim 
disk, all the gravitational energy release is not converted to the thermal radiation, but significant 
part of the energy is carried inward due to radial advection. Slim disk is geometrically thick, and 
can be much hotter than the standard disk (Kato, Fukue and Mineshige 1998; Watarai et al. 2000). 

One of the most significant observational characteristics of the slim disk is that the disk lumi- 
nosity can exceed the Eddington limit by up to ~ 10 times (Abramowicz et al. 1988; Szuszkiewicz, 
Malkan and Abramowicz 1996; Kato, Fukue and Mineshige 1998). This may be qualitatively 
understood as follows: At any disk radius, local radiation pressure may not exceed the vertical 
gravitational force, such that 

. . cGM h 

F(r) < r - 4 

kt r 

where F(r) is the energy flux, r the disk radius and h the half-thickness. Therefore, 

[Tout 

L disk = 2 2irrF(r)dr 

„ AircGM r° ut h , 

< / -?dr 

k J r . r z 
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where ri n and r ou t are inner and outer disk radius, respectively, such that In {r utl T in) ~ 10. In the 
standard disk, which is geometrically thin, h/r < 0.1, thereby Ldisk % LEdd- On the other hand, 
with the slim disk in which h/r « 1, Ldisk i$ 10 LEdd- 

In addition to that super-Eddington luminosities are permissible, slim disk has following ob- 
servational characteristics (e.g., Kato, Fukue and Mineshige 1998; Watarai et al. 2000; Watarai, 
Mizuno and Mineshige 2001): (1) As the mass accretion rate exceeds the critical rate Mq, energy 
conversion efficiency decreases due to advection (appendix A, figure 7). Hence the disk luminosity 
is no longer proportional to the mass accretion rate, but saturates at ~10 LEdd- (2) Innermost disk 
radius can be closer to the black hole than the last stable orbit, as mass accretion rate increases. 
Even in the Schwarzschild metric, innermost disk region inside 6 r g can emit significant amount of 
X-rays. Thus slim disk can be "hotter" than the standard disk. (3) If radial dependence of the disk 
effective temperature is written as T e ff oc r~ p ("p-free disk" below), the exponent p reduces from 
0.75 (which is expected for the standard disk) to 0.5 as advection progressively dominates. Thus 
spectral difference from the standard disk will be noticeable. 

4.2. Slim disk for Galactic Sources? 

Watarai et al. (2000) suggested that the slim disk model may explain the apparently hard 
energy spectra of Galactic superluminal jet sources. In fact, energy spectra of GRO J1655-40 
(Kubota 2001), GRS 1915+105 (Yamaoka 2001) and XTE J1550-56 (Kubota 2001) are not fitted 
with a standard accretion disk model when their disk luminosities and temperatures reached max- 
ima, but better fitted with the p-free disk model with p between 0.75 and 0.5. This is considered 
to be an evidence of emergence of the slim disk, when mass accretion rates are extremely high in 
these sources. 

On the other hand, when the disk luminosity and temperature are lower, GRO J1655-40 is 
considered to embrace the standard disk, since the standard disk model can fit the observed spectra 
well, and the inner disk radius is fairly constant over a large luminosity variation (Sobczak et al. 
1999; Kubota, Makishima and Ebisawa 2001). Without advection, the standard Kerr disk with 
high inclination is successful to explain the observed hard disk spectra (section 3.1). Therefore, 
except when the disk luminosity is near the peak and the slim disk presumably takes place, the 
standard disk is expected to be present in the Galactic superluminal jet sources. 

4.3. Application to IC342 Source 1 

As for ULXs, the slim disk model seems to be more likely, as ULXs are intrinsically bright 
systems. Mizuno, Kubota and Makishima (2001) studied spectral variations of several ULXs, and 
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commonly found anti-correlation between the MCD parameters Ri n and Tj„, which is another 
observational feature of the slim disk (Watarai et al. 2000). Watarai, Mizuno and Mineshige (2001) 
calculated X-ray energy spectra of slim disks, and suggested that ULXs are slim disks around 10 
- 30 Mq black holes shining more luminously than Eddington limits. A characteristic spectral 
transition has been observed in IC342 Source 1, such that the bright state in 1993 was represented 
with the MCD model spectrum, while in 2000 the source was dimmer and the energy spectrum 
was power-law like (Kubota et al. 2001). This spectral transition may be interpreted as the disk 
state transition between the blight slim disk state and another anomalous state (Kubota, Done and 
Makishima 2002). 

In the present paper, we fit the observed energy spectra of IC 342 Source 1 in 1993 with the 
same slim disk model calculated by Watarai, Mizuno and Mineshige (2001). Namely, a = 0.01, 
and the local emission is a diluted blackbody with T co i/T e ff = 1.7. The model assumes the face-on 
geometry, and the distance is 4 Mpc. The energy spectra are calculated for each grid point of M 
= 10, 32, 100 M Q and M/(L Edd /c 2 ) = 1, 3, 10, 32, 100, 320, 1000. We fitted the observed spectra 
using XSPEC spectral fitting package with this grid-model, so that XSPEC interpolates the model 
spectra for other M and M values. 

Fitting result for IC342 Source 1 ASCA GIS spectrum is shown in table 3. We find M = 
19.8M and M/Mc ~ 20. This mass accretion rate would give 20 times the Eddington luminosity 
in the case of the standard disk. However, since slim disk is less efficient at M > Mq (appendix 
A and figure 7), the disk luminosity is in fact only ~ 6 times the Eddington luminosity, which is 
allowed in the slim disk (section 4.2). 

The slim disk model can fit the observed spectra reasonably well (x 2 /dof = 1.40), but not as 
good as the standard disk models (x 2 /dof pa 0.90; table 2) (figure 4). Also, the hydrogen column 
density is significantly larger with the slim disk model (~ 9 x 10 21 cm" 2 ) compared to those with 
the standard disk models (3-5 xlO 21 cm" 2 ). This is because the spectral difference between the 
best-fit slim disk model and standard disk spectra is largest in the lower energy band (figure 5). 
We simply assumed diluted blackbody local spectra with constant T co i/T e ff, which is known to be 
a good approximation for the standard disk (section 2.3). However, precise theoretical calculation 
of the local spectrum in the slim does not exist yet, and more realistic slim disk spectral model is 
expected to fit the observed spectra better (section 4.5). 

4.4. Spectral Variation of IC342 Source 1 

The characteristic anti-correlation between Ri n and Tj„ discovered by Mizuno, Kubota and 
Makishima (2001) is considered to be a characteristic of the slim disk. Watarai, Mizuno and 
Mineshige (2001), through indirect comparison of the MCD parameters and slim disk spectra, 
claimed that this spectral variation may be explained with a slim disk at a constant M only by 
varying the mass accretion rate. We demonstrate this claim by directly fitting the ASCA IC 342 
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Source 1 spectra with the slim disk model. 

We use the same spectral datasets used by Mizuno, Kubota and Makishima (2001). The 
observation period is split into five periods depending on the flux levels. For each period, GIS and 
SIS spectra are fitted simultaneously to achieve the better statistics. In the top panel of figure 
6, we show M and M variation obtained with the face-on GRAD model. We fixed the hydrogen 
column density at the average value (Nh = 5.2 x 10 21 cm -2 ), so that the free parameters are only 
M and M. The reduced \ 2 is less than unity for all the five spectra. We see clear anti-correlation 
between these two parameters, which is just rephrasing the Ri n -Ti n anti-correlation discovered by 
Mizuno, Kubota and Makishima (2001). Since M must not vary in the real world, this relation is 
telling that the standard accretion disk is not a proper model for IC 342 Source 1. The spectral 
variation is clearly seen using only the 2-10 keV band, though less clear than using the entire 0.5 
- 10 keV band. 

Since the present slim disk model does not fit the observed spectra well below ~1.5 keV (table 
3 and figure 4), we study spectral variation only using 2-10 keV (then reduced x 2 wm De ~1). For 
fair comparison, the hydrogen column density was fixed to the same average value as used in the 
GRAD model fitting. The M and M relation obtained from the slim disk fitting is shown in the 
bottom panel in figure 6. We see that M in the range of 22.5 M & to 23.9M Q is consistent with all the 
five spectra. Therefore, based on the slim disk model, we may interpret the observed characteristic 
spectral variation of IC 342 Source 1 as a result of simple mass accretion rate variations. 

4.5. Future Problems of the Slim Disk Model 

Compared to the standard disk model, X-ray spectral study from the slim disk has only a 
short history. As already stated above, the present slim disk model may not fit the observed ASCA 
spectra perfectly well, though the spectral variation is better described with the slim disk model 
than with the standard disk. There are several difficult problems to calculate the slim disk spectra 
precisely. Although the face-on geometry is assumed in the present slim disk model, observed flux 
from a slim disk is considered to be strongly inclination angle dependent. The flux drops with 
inclination angle i more rapidly than cosi, because the horizontal photon diffusion time-scale gets 
longer than the in- fall time-scale (Kato, Fukue and Mineshige 1998). Also, relativistic effects have 
yet to be taken into account, and it is pointed out that the disk around a Kerr black hole can be 
much hotter than the Schwarzschild case (Beloborodov 1998). In addition, effects of radial photon- 
trapping may not be negligible under extremely high accretion rates (Ohsuga et al. 2002). Due to 
these effects, T co i/T e ff can be much higher than ~1.7 and may have a significant radial dependence 
(Shimura and Manmoto 2003; Kawaguchi 2003). More precise spectral model calculation from the 
slim disk is anticipated overcoming these difficulties. 
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5. Origin of the Ultra-Luminous X-ray Sources 

We have shown that the present slim disk model, though still primitive, is likely to explain 
the observed super-Eddington luminosities, hard energy spectra, and spectral variations of IC342 
Source 1, in agreement with Watarai, Mizuno and Mineshige (2001). Whereas Watarai, Mizuno 
and Mineshige (2001) indirectly compared the observed spectra and the slim disk model spectra, we 
directly fitted the observed spectra with the slim disk model, and made a quantitative comparison. 
Thereby, we obtained the black hole mass ~ 20M Q , and the disk luminosity ~ 6 times the Eddington 
luminosity for IC342 Source 1. Thus, "intermediate mass black holes" are not required to explain 
the observed luminosity ~ 10 40 erg s~ 4 . 

In our Galaxy, the most massive stellar black hole is ~ 14M Q in GRS 1915+40 (Greiner, Cuby 
and McCaughrean 2001), as far as currently measured. However, more massive stellar black holes 
may well exist in other galaxies. In fact, stellar evolution theory says main-sequence stars can have 
a maximum mass ~ 60M Q (Schwarzschild and Harm 1959), and the final black hole mass could 
be theoretically as large as its progenitor beyond 40 Mq (Fryer 1999). Observationally, Grimm, 
Gilfanov and Sunyaev (2003) found that there exists such a universal luminosity function of X-ray 
binaries that is applicable to several different galaxies including Milkyway, where normalization 
for different galaxies are proportional to the star-forming rates. Galaxies which are active in star 
formation tend to have more luminous and massive compact objects, and the maximum cut-off 
luminosity is ~ 10 40 erg s -1 . We suggest that ULXs with luminosities of ~ 10 40 erg s _1 arc 
black holes having a few tens of solar mass, which reside at the brightest end of the X-ray binary 
luminosity function. Observed evidence of the association between some ULXs and star forming 
region may support our idea (e.g., Matsushita et al. 2000). 



6. Conclusion 

1. In order to solve the "too- hot disk" problem of ULXs and Galactic superluminal jet sources, 
we have carefully studied energy spectra of the standard accretion disk and slim disk. In par- 
ticular, we have calculated the extreme Kerr disk spectral model, and studied how relativistic 
effects and disk inclination angle affect the observed disk spectra. 

2. We have found that the standard Kerr disk model can successfully explain the observe hard 
spectra of GRO J1655-40, because the Kerr disk spectra become significantly harder than 
Schwarzschild ones when the disk is significantly inclined, which is exactly the case for GRO 
J1655-40. Another super-luminal jet source GRS 1915+105 is also a highly inclined system 
(Greiner, Cuby and McCaughrean 2001), so that the idea of inclined Kerr disk model may 
be plausible too. 

3. The Kerr disk spectra are not significantly harder than the Schwarzschild disk when disk incli- 
nation angle is not large. We conclude that the standard Kerr disk model is not appropriate 



-16- 



to explain the observed hard spectra of most ULXs, since (1) it is unlikely that accretion 
disks in most ULXs are preferentially inclined, and (2) if the inclined Kerr disk is applied, 
unreasonably large black hole mass (~ 300 M ) is required. 

4. We have calculated the slim disk spectra, and applied to the observed spectra of IC342 Source 
1 in 1993, when the source is in the high state. Although the current primitive slim disk 
model does not perfectly fit the data, we found the slim disk can explain the observed super- 
Eddington luminosity, hard X-ray spectra, and spectral variation successfully. In particular, 
the observed characteristic spectral variation is explained with a constant mass at M ~ 2OM , 
only varying the mass accretion rate. 

5. We propose ULXs are binary systems with a few tens of solar mass black holes, which reside 
at the bright end of the X-ray binary luminosity function. Such moderately massive stellar 
black holes may not exist in our Galaxy, but presumably not uncommon in the galaxies where 
massive star formation is much more active. 



We are grateful to Prof. K. Makishima and Dr. Luigi Foschini for useful comments and discus- 
sion, and to Drs. D. Bhattacharya, S. Bhattacharya and A. V. Thampan for pointing out bugs in 
the old GRAD code and comparing the fixed GRAD code with their accretion disk model. We are 
glad to Dr. Hanawa for developing the original GRAD code and fixing the bugs. We acknowledge 
Dr. A. Laor for making his transfer function calculation available with the XSPEC package. We 
are thankful to the anonymous referee for many variable comments to improve the paper. This 
research has made use of public data and software obtained from the High Energy Astrophysics 
Science Archive Research Center (HEASARC), provided by NASA's Goddard Space Flight Center. 
This work was supported in part by Polish KBN grant PBZ-KBN-054/P03/2001. 



A. Energy Conversion Efficiency and Critical Mass Accretion Rates 

Energy conversion efficiency of optically thick accretion disk (rj) is dependent on disk model 
assumptions, so that a particular care is needed to interpret the disk parameters we obtain from 
spectral model fitting. In particular, conversion efficiency of the slim disk model is dependent on 
the mass and mass accretion rates (Watarai et al. 2000; Watarai, Mizuno and Mineshige 2001). 

In the Newtonian case, energy loss in the standard disk per second (= disk luminosity, Ldisk) 
may be written as 

■ / GM L 2 \ 
Ldisk = M — 2 

\ fin ^ r in J 

GMM 
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where rj n is the innermost radius, and L is the specific angular momentum at rj n , \[GMr^. If we 
take ri n as the last stable orbit in the Schwarzschild metric, 6r g = 6 GM/c 2 , Ldisk = t^^c 2 > hence 
T] = 1/12 for the Newtonian disk. In the Pseudo-Newtonian potential, which the present slim disk 
model adopts, 

■ / GM L 2 \ 

Ldisk = M ( r . _2r " 2^ j ' 

\ ' in 41 g 4i ln J 

where L = yj 'GMr^ l{r% n — 2r s ) 2 . With r{ n = 6r g , we obtain 77=1/16. This is close enough to the 
correct efficiency in the Schwarzschild metric, r/=0.057. For the extreme Kerr case with a = 0.998 
and ri n = 1.24 r g , r\ = 0.366. 

We may define the critical accretion rate Mq so that LEdd = V^cc 2 , where Lem is the 
Eddington luminosity. For the Schwarzschild case with rj n = 6r g , M c = 2.9 x 1O 18 (M/M ) g s _1 , 
and for the Newtonian case with the same innermost radius, Mc = 2.0 x 10 18 (M/Mq) g s -1 . For 
the extreme Kerr case M c = 4.6 x 1O 17 (M/M ) g s -1 . 

In the pseudo-Newtonian potential, if the innermost radius is constant at rj„ = 6 r g , Mc = 
2.7 x 1O 18 (M/M ) g s -1 . However, in the present slim disk model, the inner radius becomes smaller 
than 6r g , and the energy advection is dominant at high accretion rates (Watarai et al. 2000). Hence, 
the conversion efficiency is dependent on the mass and mass accretion rate. In figure 7, we show 
conversion efficiency of the slim disk for several mass and mass accretion rate values. We see that 
at low mass accretion rate limits, r/ ~ 1/16, but the efficiency significantly decreases with mass 
accretion rates when M > Mc due to the advection effect. 



B. Bugs in the original GRAD model 

Three mistakes have been found in the appendix of Ebisawa, Mitsuda and Hanawa (1991; 
EMH) besides obvious typos: 



(A8) in EMH should be read as follows: 

d0' ph cos i sin 6* p h + sin i cos # p h sin ^ p h 



d8- 



P h 



sin#' 



P h 



(A15) in EMH should be read as follows: 



(Bl) 



3GM 



rc* 



1 



2GM_y 05 /GMrdip ph 
rc 2 J V rc 2 c dt 



P h 



The bold-face part was 1 in EMH, which is wrong. 
(A22) in EMH should be read as follows: 



T, 



col 



1.11 



T co i/T c 



off 



1.5 



g(r/r g ) 
0.15 



M 



1.0 M Q 



-1/2 



M 



1/4 



(B2) 



10 18 g s" 



keV. 



(B3) 
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The bold-face part was 1.4 in EMH, which is incorrect. 

Below, equations (Al) to (A4) refer to those in the present appendix, not in EMH. Equation 
(Bl) was already fixed in the GRAD model code used in EMH, which was distributed with XSPEC 
version 11.0. lae and previous ones. Equations (B2) and (B3) were not fixed in the code, and the 
wrong values were used in EMH. In addition, there was another bug in the code used in EMH: In 
the part to calculate the disk flux, 

F = 8.9038633D - 3 * F/Ratio * *4/(D/10.0) * *2 * (Em/1.4) * *2 

was wrong, and this has to be 

F = 8.9038633D - 3 * F/Ratio * *4/(D/10.0) * *2 * Em * *2, (B4) 

where Em denotes the compact object mass. These bugs in equation (B3) and (B4) compensate 
each other (but not perfectly), which probably explains why we could not notice these bugs much 
earlier. The bug in equation (B2) affects only when the disk is inclined, and this effect is not 
significant. 

Besides the bug in equation (B2), one can see from equations (B3) and (B4) that the old 
GRAD spectra with M = 1.4M and the new GRAD spectra with M = l.OM are identical. 
Therefore, besides the minor difference of redshifts which originates in (B2), when an observed 
spectrum was fitted with the old GRAD model, the mass of the compact object was erroneously 
1.4 times higher than the correct mass given by the fixed GRAD model. This should be taken into 
account to interpret published results obtained using the old GRAD model (e.g., Ebisawa, Mitsuda 
and Hanawa 1991; Kitamoto, Tsunemi and Roussel-Dupre 1992; Ebisawa et al. 1993; Dotani et al. 
1997; Kubota et al. 1998). The mass accretion rate represented in the dimension of mass per time 
is unchanged. The fixed GRAD code is distributed with XSPEC version 11.1.0 and on. 

Note that the bug in the old GRAD code may compensate with an ambiguity of T co i/T e ff 
values. In fact, in Ebisawa, Mitsuda and Hanawa (1991) and Ebisawa et al. (1993), we used 
T co i/T e ff = 1.5 and obtained reasonable masses of black holes and neutron stars, which should 
be reduced by 1.4 times using the correct GRAD model. On the other hand, the mass estimate 
is proportional to (T co i/T e ff) 2 (when distance is known) or (Tcoi/Teff) 4 " (when d/y/M is known), 
thus using T co i/T e ff = 1.8 or 1.6 instead of 1.5 retrieves the original mass values. The higher 
Tcoi/Teff values are in fact consistent with a precise calculation carried out more recently (Shimura 
and Takahara 1995; Ross and Fabian 1996). 

We confirmed that the fixed GRAD model gives almost identical spectra with the Schwarzschild 
disk spectra calculated by Bhattacharyya, Bhattacharya and Thampan (2001). Only a minor dif- 
ference is found when the accretion disk is more inclined than ~ 60°, when the model by Bhat- 
tacharyya, Bhattacharya and Thampan (2001) gives as much as ~ 10 % higher flux at ~ 10 keV. 
The discrepancy is probably due to the fact that the GRAD model does not take account of the 
photons whose trajectories change directions more than 180 degrees. In fact, those photons may 
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not be observable being blocked by the putative plasma near the black hole which is responsible 
for generation of the power-law hard-tail component. 



C. Comparison of Newtonian and Schwarzschild Accretion Disk Spectra 



Comparison of the GRAD and MCD model parameters have been made by EMH, which is 
affected by the bug explained above. Here, we present correct results using the new GRAD code 
after the bug is fixed. Also, we discuss how the relativistic effects (in the Schwarzschild case) affect 
the mass and mass accretion rate determination from observed accretion disk spectra. 

MCD model does not take into account the inner boundary condition, but simply assumes 
T(r) = T in {R in /rf/ A , where R in and Ti n are the apparent inner disk radius and temperature, 
respectively. MCD model is useful to represent observed spectra, since it has the two independent 
parameters which are directly constrained from observation, such that the spectral shape is deter- 
mined only by Tj n , and the normalization is proportional to R^ n cosi, where i is the inclination. 
Often disk inclination is unknown, so that Ri n is uncertain and proportional to (cosi) -0 ' 5 . 

Relationship between the MCD Tj n and R{ n parameters and the mass and mass accretion rate 
in the Newtonian disk model with correct boundary condition is given as (see also Kubota et al. 
1998), 




(Ri n km) M Q , 
1/4 'MN — 



Mr. 



© 



keV. 



(CI) 
(C2) 



The multicolor disk spectrum and the Newtonian disk with correct boundary condition have 
almost identical spectral shape in 0.5 - 15 keV when equations (CI) and (C2) are considered. 

Formulae to relate the new GRAD model parameters and MCD parameters are the following: 
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The numerical coefficients have been obtained by fitting the new GRAD spectra with multicolor 
disk model (see below) 3 . 

In Figure 8, we illustrate the relationship of the MCD model parameters and M and M in the 
Newtonian and GRAD models. These figures are made by simulating GRAD and Newtonian model 
spectra with various values of mass, mass accretion rate and inclination, using the Ginga response 
(sensitive in 1 - 30 keV), assuming the distance 1 kpc and exposure time 3000 sec exposure. These 
simulated spectra are fitted with the MCD model, and the best-fit Ri n and Tj„ are obtained. 

From equations (CI) to (C4) and Figure 8, for a given observed spectrum for which Ri n and 
Ti n are determined, we can see the following (see also table 1 and 2): 

• When fitting the same observed spectra, we will get a smaller mass and higher accretion rate 
with the GRAD model than with the Newtonian disk. 

• Differences in the mass and mass accretion rate between the GRAD model and the Newtonian 
models are the largest for the face-on disk, and get smaller with inclination. 

These characteristics of the GRAD model are due to relativistic effects, and may be understood 
intuitively as follows. Because the conversion efficiency of the GRAD model (0.057) is smaller than 
that of the Newtonian disk (0.083), the GRAD model requires a higher mass accretion rate to shine 
with the same luminosity. Thus we need a higher mass accretion rate with the GRAD model when 
fitting the same observed spectra. The observed disk flux is proportional to L^isk cos i, thus the disk 
luminosity induced from the observed flux is L^isk °c (cosi) -1 . Consequently, the mass accretion 
rate is proportional to (cosi) -1 in the Newtonian case. In the GRAD model, this dependence 
becomes milder, because the face-on disk tends to be dimmer due to the light bending effect in the 
vicinity of the black hole, while those photons removed from the near face-on disks turn out to be 
the flux enhancement for the inclined disks. In fact, for the face-on disk, the GRAD mass accretion 
rate is ~ 1.8 time larger than that of the Newtonian disk, whereas it is ~ 1.1 times for i = 80°. 

Because of the gravitational redshift, the GRAD model gives a smaller characteristic disk tem- 
perature than the Newtonian disk with the same mass and mass accretion rate (compare equations 
C2 and C4). To achieve the same characteristic disk temperature of the Newtonian model, the 
~1.8 times mass accretion rate increase is not sufficient, and the mass has to get smaller too. Since 
we measure the projected disk area oc M 2 cosi, the mass increases with inclination as (cosi)" ' 5 in 
the Newtonian case, while the GRAD disk temperature gets higher because of the Doppler boost. 
Therefore, ratio of the GRAD mass to the Newtonian mass is the smallest for the face-on disk 
(~0.6) and gets larger with inclination (~0.73 for i > 60°). 

It is interesting that difference of the Newtonian and GRAD model parameters are the largest 



'^Equations (3) in Ebisawa, Mitsuda and Hanawa (1991) was wrong such that the old formula gives 1.4 times larger 
mass for the same Ri n and Ti„. 
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for the face-on disk and the smallest for highly inclined disks where relativistic effects are supposed 
to be most significant. This is because the gravitational redshift and Doppler boosts more or less 
cancel each other for a highly inclined Schwarzschild disk. This is in contrast to the extreme Kerr 
disk, where near edge-on disks show extremely hard spectra due to enormous Doppler boosts (see 
figure 1 and 2). 
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Table 1. GRO J 1655-40 fit with Newtonian, Schwarzschild and Kerr disk models' 



Newtonian disk Schwarzschild disk Kerr disk (a — (1,!)<)S) 
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M/Mq 


Af b 


M/M c ° 


N pom d 


N H ° 


x2/dof 


M/M e 


Af b 
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N d 
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N pom d 




x2/dof 


0° 


1.3 


0.49 


0.18 


7.2 


0.97 


1.3 


0.83 


0.87 


0.37 


6.9 


0.97 


1.35 


2.6 


0.58 


0.48 


6.2 


0.98 


1.40 


30° 


1.4 


0.56 


0.19 


7.2 


0.97 


1.3 


0.96 


0.99 


0.36 


6.0 


0.95 


1.25 


3.0 


0.33 


0.24 


4.0 


0.94 


1.20 


45° 


1.6 


0.69 


0.22 


7.2 


0.97 


1.3 


1.2 


1.2 


0.34 


5.2 


0.93 


1.17 


5.2 


0.36 


0.15 


3.9 


0.94 


1.30 


70° 


2.3 


1.42 


0.31 


7.2 


0.97 


1.3 


1.8 


2.1 


0.40 


4.7 


0.92 


1.08 


15.9 


0.35 


0.048 


8.9 


0.99 


1.45 


80° 


3.2 


2.80 


0.43 


7.2 


0.97 


1.3 


2.4 


3.1 


0.45 


6.3 


0.93 


1.14 


27.0 


0.39 


0.032 


10.4 


0.99 


1.62 



a Distance is assumed to be 3.2 kpc, and T coi /T e j.- j — 1.7. Slope of the power-law component is fixed to 2.5. One percent systematic error is 
included for each spectral bin. Degree of freedom (dof) is 156. 

b Unit is 10 1S g s _1 . 

c The critical mass accretion rate is so denned that M = M<j gives the Eddington luminosity. See Appendix A. 

Power-law component normalization in photons s cm ' 2 kcV -1 at 1 kcV. 
°In 10 22 cm" 2 . 
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to 
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Table 2. IC342 Source 1 fit with Newtonian, Schwarzschild and Kerr disk models a 



Newtonian disk Schwarzschild disk Kerr disk (a — 0.998) 
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M/Mq 


M b 


M/M c c 
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M/Mq 
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X 2/dof 


M/Mq 
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0.50 
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0.42 


0.91 


10.1 


335 
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0.44 


0.90 


29.5 
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0.51 


0.90 
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0.41 
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387 


11.2 


0.44 


0.90 


51.3 


113 


4.8 


0.51 


0.90 


60° 
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341 
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0.42 


0.91 


15.2 


502 


11.4 


0.44 


0.89 


107 


120 


2.4 


0.46 


0.91 


80° 


35.1 


981 


14.0 


0.42 


0.91 


25.9 


1069 


14.2 


0.38 


0.91 


332 


153 


1.0 


0.30 


1.00 



a Distance is assumed to be 4 Mpc, and T co i/T e ff — 1.7. Degree of freedom (dof) is 52. 
b Unit is 10 18 g s _1 . 

c The critical mass accretion rate Mc is so defined that M = Mc gives the Eddington luminosity. See Appendix 

A. 

d In 10 22 cm" 2 . 



Table 3. IC342 Source 1 fit with Slim disk model a 



i M/Mq 


M b 


M/M c c 


Ldisk/LEdd 


N H d 




0° 19.8 ±1.7 


iooq±S88° 


20±f3 


5.9 


0.90±8;8I 


1.40 



a Distance is assumed to be 4 Mpc, and T co i/T e ff =1.7. The viscous parameter a = 0.01. Degree 
of freedom (dof) is 52. Errors correspond to 90 % confidence level. 

b Unit is 10 18 g s _1 . 

c The critical mass accretion rate Mq is so defined that M = M<j gives the Eddington luminosity 
in the case of the standard disk. See Appendix A. 

d In 10 22 cm" 2 . 
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Comparison of Newtonian, Schwarzschild and Kerr disk spectra 
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Fig. 1. — Comparison of Newtonian (black), Schwarzschild (red) and extreme Kerr (a = 0.998; 
green) disk spectra for the face-on (i = 0°) and a near edge-on disk (i = 80°). Inner disk radius 
is 6 r g for the Newtonian and Schwarzschild disks, and 1.24 r g for the Kerr disk. The same mass 
accretion rate is assumed, which is so chosen to give the Eddington luminosity for the Schwarzschild 
disk. Other disk parameters are indicated in the figure. Note that the total disk luminosities are 
different for the three disk models because the energy conversion efficiencies are different (0.083, 
0.057 and 0.366 for Newtonian, Schwarzschild and Kerr disks, respectively). 
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Kerr disk at d = 1 kpc, M= 1 M Q , Eddington limit 
/x = 0.9, 0.7, 0.5, 0.3 and 0.1 

o 

o 

o b I 1 1 1 1 1 1 1 1 1 




o 0.01 0.1 1 10 



Energy [keV] 

Fig. 2. — The Kerr accretion disk spectra with an extreme angular momentum (a = 0.998), 
observed at the inclination angle fj> = cosi = 0.9 (green; near face-on), 0.7 (yellow), 0.5 (cyan), 
0.3 (red) and 0.1 (black; near edge-on). Note the unit of the ordinate (keV 2 s" 1 keV -1 cm" 2 ) 
which facilitates to see the energy release per logarithmic energy. Solid lines indicate the total disk 
spectra, and contributions from inner (1.26r g < r < 7r g ), middle (7r g < r < 400r 9 ), and outer 
parts (400r 9 < r) are plotted separately either by dotted line or broken line. The distance and 
mass are assumed to be 1 kpc and 1 M Q respectively. The Eddington luminosity is assumed, and 
T co i/T e ff = 1. 
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GRO J1655-40 IC342 Source 1 



T 




Energy [keV] Energy [keV] 

Fig. 3.— ASCA GIS energy spectra of GRO J1655-40 (left panel) and IC342 Source 1 (right). 
The best-fit Schwarzschild disk models are indicated in red, in which case the best-fit masses 
are too small (1.8 M Q and 8.9M Q for GRO J1655-40 and IC342 respectively) compared to the 
mass determined from optical observations (7 Mq for GRO J1655-40), or that expected from the 
observed luminosity (100 Mq for IC342). Schwarzschild disk spectra with the expected masses 
give too low temperatures to explain the observed accretion disk spectra for both cases (shown in 
green). Additional power-law component (blue in left panel) is required to fit the GRO J1655-40 
spectrum. 



-31- 



IC 342 Source 1; Standard Disk Fit 




Energy [keV] 

Fig. 4. — IC342 Source 1 ASCA GIS spectrum fitted with the Schwarzschild disk model (i = 0°; 
top) and with the slim disk model (bottom). Residual is more conspicuous in the slim disk model 
below ~ 1.5 keV. 
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Fig. 5. — IC342 Source 1 best-fit standard accretion disk model (face-on Schwarzschild disk; black) 
and slim disk model (red). Models removed of interstellar absorption are shown in dotted lines. 
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Standard disk model 
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Fig. 6. — Spectral variation of IC342 Source 1 observe in September 1993. The same datasets 
as in Mizuno, Kubota and Makishima (2001) are used. The upper-panel indicates the variation 
of mass and mass accretion rates obtained by applying the standard Schwarzschild disk model 
(GRAD model). Contours with dotted lines are from fitting in 0.5-10 keV, and those with solid 
lines are from 2-10 keV. The two contour levels indicate 68 % (1 a) and 90 % error regions for 
two parameters. The lower panel is obtained through application of the slim disk model in 2 - 10 
keV. The region marked with gray indicates that the observed spectral variation is achieved with 
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Fig. 7. — Energy conversion efficiency of the slim disk used in the present paper and in Watarai, 
Mizuno and Mineshige (2001). The vertical broken lines indicate the critical mass accretion rates 
Mc which gives the Eddington luminosity in the standard disk model (see appendix A). 
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Fig. 8. — Comparison of the Newtonian and Schwarzschild disk spectra through representing these 
spectra with MCD model parameters. Upper figure shows Ri n and Tj n values to fit each of the 
Newtonian or Schwarzschild disk spectra with various mass, mass accretion rate and inclination. 
Note that Ri n is proportional to mass. In the lower figure, ordinate is R^ n T^ n , which is proportional 
to the disk luminosity and mass accretion rate. Unit of the mass accretion rate is 10 18 g s _1 . 
Distance is assumed to be 1 kpc, and T co i/T e ff = 1. 



